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Recently, generalizations of fractional quantum Hall (FQH) states known as fractional quantum 
anomalous Hall or, equivalently, fractional Chern insulators states have been realized in lattice 
models. Ideal wavefunctions such as the Laughlin wavefunction, as well as the corresponding trial 
Hamiltonians for which the former are exact ground states, have been vital to characterizing FQH 
phases. The Wannier function representation of fractional Chern insulators proposed in [X.-L. Qi, 
Phys. Rev. Lett. 107, 126803] defines an approach to generalize these concepts to fractional Chern 
insulators. In this paper, we apply the Wannier function representation to develop a systematic 
pseudopotential formalism for fractional Chern insulators. The family of pseudopotential Hamilto- 
nians is defined as the set of projectors onto asymptotic relative angular momentum components 
which forms an orthogonal basis of two-body Hamiltonians with magnetic translation symmetry. 
This approach serves both as an expansion tool for interactions and as a definition of positive 
semidefinite Hamiltonians for which the ideal fractional Chern insulator wavefunctions are exact 
nuUspace modes. We compare the short-range two-body pseudopotential expansion of various frac- 
tional Chern insulator models at filling /i = 1/3 in phase regimes where a Laughlin-type ground state 
is expected to be realized. We also discuss the effect of inhomogeneous Berry curvature which leads 
to components of the Hamiltonian that cannot be expanded into pseudopotentials, and elaborate on 
their role in determining low energy theories for fractional Chern insulators. Finally, we generalize 
our Chern pseudopotential approach to interactions involving more than two bodies with the goal 
of facilitating the identification of non-Abelian fractional Chern insulators. 



I. INTRODUCTION 

The field of topological insulators (TIs) is currently 
witnessing enormous interest in condensed matteiQ"^!, 
The predecessor of TIs is the integer quantum Hall effect 
(IQHEjP realized in a two-dimensional electron gas with a 
strong perpendicular magnetic field. Similar to other TIs 
discovered more recently, the IQHE is a gapped state of 
matter characterized by topologically robust edge states 
and a bulk topological invariant known as the Chern 
number or the Thouless-Kohmoto-Nightingale-den Nijs 
(TKNN) numbeJ'^. The IQH state has been generalized 
to lattice models without orbital magnetic fielcP, which 
are named as quantum anomalous Hall (QAH) states or 
Chern insulators (CI). QAH states have been proposed 
in realistic materiala^lSl. In 1983, the fractional quantum 
Hall effect (FQHE) was dis covere d in systems with a frac- 
tionally filled Landau levePE^. Since then, the FQHE 
has become a. paradigmatic example of a topologically 
ordered phase^ that exhibits a non-analytic role of in- 
teractions. Numerous fundamental developments evolved 
out of this direction, such as the concept of non-Abelian 
statistics which f orms t he foundation of topological quan- 
tum computatiorP^ini. As different to IQHE fundamen- 
tally, the flatness of the Landau level and its associated 
freezing of kinetic energy appears to be necessary condi- 
tions for the FQHE state to be energetically preferable. 

This immediately provokes the question of whether a 
fractional Chern insulator (FCI) can be realized, i.e. a 
lattice version of the FQHE without an external mag- 
netic field. Taking a CI such as Haldane's honeycomb 
modei^ as a natural starting point, the task then is 
to drive the system into the fiat band limit where the 



chemical potential lies within this band, e.g. at frac- 
tional one third filling, which is well separated from the 
other bands and hence accomplishes a FQHE-type lat- 
tice scenario. (Unlike the quantum Hall case, the FCI 
filling is not given by the ratio of electrons over mag- 
netic fiux quanta, but the chemical potential of the lat- 
tice model.) Different groups have recently indepen- 
dently pursued this direction, proposing FCI models on 
the honeycomb, kagome, square, and checkerboard lat- 
ticepHUl. In different ways, the flattening of the Chern 
band can be accomplishe d thro ugh geometric frustration 
(e.g. long-range hoppingl ^^ l ^^ l, multi-band effect^, and 
multi-orbital characteiflSl. While previous s and p-type 
orbitals involved in topological insulator candidate ma- 
terials would assume only moderate interactions from 
small hybridizations, d-orbital-type systems provide an 
arena for both strong correlations and topological band 
structureJi^. First numerical investigations of the FCI 
phases on a torus at one third band filling found indi- 
cations of a three-fold topologically degenerate ground 
state separated from the other energy levels by a gap, 
where the fiux insertion showed level crossings with no 
level repulsion between them, and the Chern numb ers of 
these many-body ground states found to be 1 /3 eactP^'^ol 
While this already gives a strong hint that a Laughlin- 
type fractional Chern phase might be realized, this does 
not yet completely rule out a competing charge density 
wave (CDW) state at this filling, which can show simi- 
lar fractional Chern numbers in the ground states, level 
degeneracy, and a gap. But further evidence against a 
CDW state has been found by finite size scaling, entan- 
glement measures, and the distribution of ground state 
momenta as a function of cluster siz^^lJ, Compared to 
the FQHE for which the joint perspective of energy and 
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entanglement measures generally gives a consistent and 
complementary picture, the current stage of FCI models 
particularly calls for further investigation. 

In this article, we focus on further developing the un- 
derstanding of FCI phases from the perspective of en- 
ergetics and interactions. In general, FCIs involve dif- 
ferent scales such as the kinetic bandwidth of the frac- 
tionally filled Chern band, the gap separation from other 
bands, as well as the magnitude and range of interactions. 
Even if we assume a conventional FQHE-type parameter 
window where inter-band scattering is neglected and the 
bandwidth of the fractional Chern band is assumed small 
versus the interaction strength, a crucial complication of 
FCI models is the inhomogeneous Berry curvature which 
has no quantum Hall analogue. As a signature of this 
difference, the Platzman-Girvin-MacDonald algebrgP^I of 
the lowest Landau level (LLL) can only be mapped to 
the lattice Chern density operators in the continuum 
limit as well as for homogeneous Berry curvatur^^, from 
where Hamiltonian theories can be constructecP^. How- 
ever, an exact one-to-one mapping between FQH and FCI 
states with Chern number C = 1 has been established by 
the Wannier state representation of Chern insulatorg^S 
despite of the inhomogeneous Berry curvature. (One- 
dimensional) Wannier states are single particle states 
which are localized in real space in one spatial direction 
(such as x), but are momentum eigenstates in the orthog- 
onal direction y. Different Wannier states are related by 
translation in x direction, and all Wannier states form a 
complete basis of the single particle Hilbert space defined 
by a non-degenerate energy band. The exact mapping 
between FQH and FCI is obtained by mapping Landau 
level wavefunctions in the Landau gauge to the Wan- 
nier states in FCI. More details of this mapping will be 
reviewed in Sec. |ll] From the Wannier state represen- 
tation, we learn that the effect of inhomogeneous Berry 
curvature is absent if we consider a special Hamiltonian 
obtained by mapping a FQH Hamiltonian to the FCI sys- 
tem. In other words, the effect of inhomogeneous Berry 
curvature in an FCI strongly depends on the interaction 
Hamiltonian. In two recent worka^MIi^ the Wannier state 
representation and its further improvement has been in- 
vestigated in both fermionic and bosonic FCI systems. In 
the bosonic z/ = 1/2 FCI, the wavefunction proposed by 
the Wannier state representation has a high overlap with 
the exact ground state wavefunction obtained by exact 
diagonalization^^. In the fermionic u — l/i case, such 
a high overlap is also achieved given that the Wannier 
states are modifiecPSI. Therefore the validity of the Wan- 
nier state representation has been demonstrated at least 
in those simplest FCI states. 

In this article, we develop a pseudopotential (PP) for- 
malism for fractional Chern insulators. Previously, PPs 
have been established in the context of FQHEpSl. ppg 
are partial wave expansions of the Coulomb interactions. 
The resulting expansion quantum number is the relative 
angular momentum m of two particles (with m even for 
bosons and odd for fermions) where the expansion co- 



efficients Vm denote the energy penalty of two particles 
having a relative angular momentum of m. In the same 
way the Landau level wave functions were used as a basis 
for defining such an expansion in FQHE, we now employ 
a similar construction for the FCI Wannier functions^S. 
PPs proved extremely useful in FQHE not only to give 
a universal classification of different interaction profiles, 
but also to obtain an adequate description of general 
FQHE phase diagrams. Furthermore, many paradig- 
matic FQH wave functions are exact ground states of cer- 
tain PP Hamiltonians for which representative finite size 
studies would be more accurate to resolve their universal 
properties than for a generic interaction scenaricPsHUl. In 
adapting this concept to FCI models, we can hope for a 
similarly promising route to a deeper and more universal 
understanding of FCI Hamiltonians. 

The paper is organized as follows. In Section |llj 
we review and expand the description of the Wannier 
state representation of fractional Chern bands. This lays 
the foundation for the definition of FCI PPs, for which 
we first review the PP formalism for FQHE and sub- 
sequently develop the FCI formulation of PPs in Sec- 
tion |III[ In Section [TV] we apply our PP formalism to dif- 
ferent FCI models and their interaction profiles. We find 
that, as a direct consequence of inhomogeneous Berry 
curvature, there is a portion of the interactions which 
cannot be expanded in PPs. These observations will be 
analysed in detail in Section [V] There, we classify what 
type of center-of-mass (CM) breaking and magnetic- 
translation-group (MT) breaking scattering elements ap- 
pear in FCI interactions, and how these symmetries are 
potentially reemergent in an effec tive low energy theory 
of the problem'^^ In Section VI we generalize the 
FCI PP principle to many-body interactions, which en- 
ables us to define exact Hamiltonians for non-Abelian 



FCI phases. In Section [VII| we conclude that the pseu- 
dopotential formalism establishes a suitable platform to 
further investigate and analyze new states of matter in 
FCIs. 



II. WANNIER STATE REPRESENTATION OF 
FRACTIONAL CHERN BANDS 

In this section, we review the Wannier state represen- 
tation of FCIs proposed in Ref. [25l The idea of this 
approach is to find a suitable single-particle basis, the 
one-dimensional (ID) Wannier state basis, and to use 
this basis to establish an exact mapping between FCI 
and FQH states. While a cylindrical geometry was em- 
ployed in Ref. HSJ the discussion can also be formulated 
on the torus geometrjEU which we choose to use in the 
following. (The torus formulation of the Wannier state 
representation has also been investigated independently 
in ReiW^.) 

Consider a band insulator with the Hamiltonian 
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with i, j being the site indices of a two-dimensional lattice 
with periodic boundary conditions and a,P = 1,2,..^N 
labeling internal states in each unit cell such as orbital 
and spin states. With translation symmetry h"^^ = 

h"^ , , the Hamiltonian can be written in momentum 
space as 



H 



(2) 



with 



I t t E "^^"^ 
^ V \^ 



-ik-Ti 



Q/3g-ik-(rj-ri 



We use Lx, Ly to denote the number of lattice sites in x 
and y direction, respectively. The momentum k takes 

values of (^x^'^X^)' '^ith = Lj., Uy — 

1,2, integers. The Hamiltonian matrix /it can be 

diagonalized to obtain the eigenstates 



/ik |n,k) = £'„k I"-, k) , n^l,2, ..,N. 



(3) 



We are interested in the system with a lowest energy band 
occupied, and a gap separating this band from all 
other bands. Since only the lowest band will be involved, 
we will denote 1 1 , k) by | k) for simplicity. 

In the thermodynamic limit L^^Ly — )■ oo, k is a 
good quantum number and the Berry's phase gauge field 
a. — —i (k| Vk |k) can be defined, which determines the 
first Chern number as the flux of the gauge field in the 
Brillouin zone: C*! = ^ J^^ d^kV x a. For the realiza- 
tion of FCIs, we are interested in a band with Ci ^ 0. 
More specifically, in this paper we will focus on Ci = 1 
systems. Moreover, for finite Lx,Ly, it is necessary to 
generalize the definition of a Berry's phase gauge field 
and Chern number to the case of |k) with a discrete k 
variable. 

We start from the definition of ID Wannier states 




(4) 



which is a Fourier transform of |k) in the x-direction, but 
which remains an eigenstate of ky. Since the state |k) is 
only determined by the Hamiltonian up to a phase, the 
phase factor e*^'' is not pre-determined. As was discussed 
in Refs. [55] and ESI the phase ambiguity can be fixed by 
defining the projected position operator x = PxP with 
P — ^jj. |k) (k|, the projection operator to the occupied 
band, and x = ^ Xi \i, a) {i, a\ the position operator. 
However, in a system with periodic boundary conditions 
in a;-direction, it will be slightly more problematic to ap- 
ply this definition of x due to the dependence on the 
choice of the boundary site. As pointed out in Ref. [Ml 




FIG. 1. (Color online) Shift of the Wannier polarization as 
ky varies. The azimuthal direction represents the real space 
X while the poloidal direction represents the periodic domain 
of ky. For any full rotation along y, i.e. ky ^ ky + 2n, the 
Wannier state shifts by a; — ^ x -|- 1. 



this problem can be resolved by defining a unitary oper- 
ator 



X — exp 



2tt 



= exp 



E2nx.i X - 



(5) 



This definition preserves the periodicity x — > x + L^. The 
eigenstates of this operator are the states localized on a 
given site n in the x-direction. We define the projected 
operator 



X = PXP. 



(6) 



In momentum space, (k| X |k') = (k| X |k'). It is easy to 
see that X shifts the momentum k^ by 2tt/Lx, since 



X|k) =^^.kae'"'('=== + *^)N,a) 



(7) 



with Ukct = \/LxLy (0, a I k) the periodic part of the 
Bloch wave function. Therefore, the only nonzero ma- 
trix element of (k| X |k') is 

Fk^ky = {K + 27r/L^, ky\ X jfe^,, ky) 

a 

In the subspace of states with given ky, the matrix of X 
in momentum representation is 







X = 



F, 



2it/L^ 



\L^-l)Tr/L^ 



F2^ \ 



/ 



(9) 



where the omitted index ky is the same for all states. 



In the thermodynamic limit — >■ oo, F/ 



1 



ij^ax(k), with Qx the x component of the Berry's phase 
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gauge field. For finite L^, \Fk^ky\ should be close to 
but not exactly equal to 1. Therefore, X is an approxi- 
mately unitary matrix. To define the maximally localized 
Wannier states, we deform the X operator to a unitary 
operator by defining 



( 



X = 







-iA 



-iA 



with 



with Pxiky) 



1 

2^ 



Art 



/ ^ - -Px fey 

0<Px<27r 



This definition is periodic in 
sponding eigenvalues are 



— !• /ct + 27r. The corre- 



X\Wnk,)^e 



large limit, 



2ir , 



J2TT[P^{2Trn/Ly)-P^{2TT{n-l)/Ly)] 



(a) 



ky/2Tr 




Ky/2W 



Ky/2T, 



\ FIG. 2. (Color online) Wannier center evohition from center 
a; = to a; = 10 as a function of (a) ky/2Ti and (b) of the 
extended wave vector Ky/2Ti. From (a) to (b), it becomes 
visible how the evolution of x changes as a function of Ky, 
suggesting its similarity to Landau level wave functions (c). 



/ 

(10) 



Here, the index of the rows and columns are = 
0, -j^, 27r — j^. The eigenstates of the X operator 
form an orthogonal complete basis. Due to the simple 
form of X in momentum space, one can prove that the 
eigenstates of X are Wannier states defined in Eq. ([4]), 
with the phase (/Jk defined by 

■^k = - X! ^P'^^y ^ ^xPx{ky), (11) 

0<Px<fem 



(12) 



Therefore, we see that the center-of-mass (CM) position 
of the state jWrifc^^) is shifted by Px away from the lat- 
tice site position n. This fact indicates that Px{ky) has 
the physical meaning of charge polarization"^^. In the 



Since Px{ky) is a C/(l) phase for each fcy, one can define 
its winding number when ky goes from to 27r: 



(13) 

The log function is defined to take the value of the phase 
difference 27r [Px{2-Kn/ Ly) — Px{2TT{n — 1)/Ly)] in the re- 
gion of [— 7r,7r). As long as Ly is not too small so that 
Pxiky) would not jump by tt between two neighboring ky 
values, the Ci obtained agrees with the Chern number in 
the large Ly limit. 

For Ci ^ 0, the Wannier states \Wniky)) have a 
"twisted boundary condition" in ky, since Px(ky -|-27r) = 
Px{ky) - Ci, such that \Wn[ky + 2t:)) = \Wn+cAkv))- 
As is illustrated in Fig. [TJ for Ci = 1 the Wannier state 
CM position Xn{ky) = n+Px{ky) forms a helical curve on 



the parameter space torus of x,ky. The key observation 
which enables the Wannier state representation of FCI 
is the fact that the twisted boundary condition allows to 
label all Wannier states in such a 2D system by a ID 
parameter. If we define 



+27™ 



)^\WnMy), iorkye[0,2Tr), (14) 



it yields |lFx), with K — ky + 2iin a continuous function 
of K (in the large Ly limit). The CM position xk — 
{Wk\X\Wk) = e*T7("+^-(*=«» for K e [2nn,2n{n + 1)) 
is continuous in K and satisfies XK+2Tr — xk -I- 1. In this 
sense, xk increases linearly with K (Fig.[2|. 

Due to this behavior of an exact mapping can 

be defined between the Wannier states in the Ci = 1 FCI 
and the LLL states in a FQH problem. Consider a spin- 
less fermion with the Hamiltonian H = 2^(p — A)^ on a 
torus of the size LxIb x Lylg with the uniform perpen- 
dicular magnetic field V x A = B = 2t: The total 
number of fiux is -/V$ = LxLy, so that the LLL contains 
the same dimension of Hilbert space as the lattice model 
discussed above on a lattice of the 



gauge choice 
have the form 

ipK{x,y) 



0, Ay = Bx, the LLL wave functions 



E 



ez 



(15) 



with 'd{z\T) the Jacobi theta functiorP^ which are appro- 
priate superpositions of the cylindric wave functions. 

Notice that we have defined the momentum K slightly 
differently from the usual definition used in Ref. [5S] and 
[551 and in the subsequent sections, so that here K is 



dimensionless and given by K 



27r „ 



n S Z on the 



LxIb X Lyls torus. This definition leads to identical re- 
sults as the usual definition used later if we replace the 
Ib here by y/2nlB- 

For Lx ^ 1, this wave function is a Gaussian function 
around the CM position xk = K^^- Denoting \4'k) as 
the state corresponding to wave function ipxix^y), one 
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can define a unitary mapping between the Hilbert spaces 
of tlie Landau level and the lattice Ci Chern insulator: 



III. PSEUDOPOTENTIAL REPRESENTATION 
OF INTERACTIONS 



/: Hci 

fi\WK)) 



Hlll, 



(16) 



with Hci and Hlll denoting the Hilbert spaces of the 
CI and LLL, respectively. Such a mapping preserves 
the continuity in K and also the topological properties 
of \Wk} and IV'if): i-e., their winding while momentum 
K is increased. Using the reverse map f~^, the many- 
body states of the LLL, such as Laughlin states and all 
other FQH states defined in the LLL, can all be mapped 
to corresponding states in the FCL Similarly, a Hamil- 
tonian iJ of a FQH system can also be mapped to a 
corresponding Hamiltonian -ffpci = f~^Hf- The main 
purpose of the current paper is to study the Hamiltoni- 
ans Hpci which are mapped from the PP Hamiltonians 
in the FQH system. One can also perform the reverse, 
mapping the FCI Hamiltonian such as a Hubbard type 
interaction Hjj of the lattice fermions to a FQH Hamilto- 
nian HpQn = fHuf~^. More details of such a mapping 
will be evaluated in the following sections. 

There is a subtle point that we want to discuss before 
ending this section. The definition of maximally local- 
ized Wannier states in Eq. [TT] still leave an ambiguity in 
the relative phase between different |W„fc^). If we rede- 
fine \Wnky) e"^"^ l^riky) with any phase Ok^, all the 
properties discussed above remain the same. The map /, 
however, depends on this choice and different choices of 
phase corresponds to physically different mappings be- 
tween FCI and FQH systems. To preserve the locality 
in the mapping, a choice should be made which makes 
jW^nfej,) continuous in ky in the large Ly limit. An exam- 
ple of the choice is the foUowin^Sl 



PyO - - 



0<Py<ky 

— V 4' 

9^ Py ' 



(17) 



Py 



0<py<2-n- 
\ a 



+ ^,a'^0,py,a 



In the L,,, 



oo limit, A' 

' Py 



ay{0,Py) 



2tt 

Ly 



This choice of 



9ky corresponds to a gauge transformation which makes 
aj,(k) uniform along the — line. Any other gauge 
choice 9'^. = 6ky + S9ky also works and describes topo- 
logically equivalent states, as long as S9k is a smooth 
periodic function of ky in the large Ly limit. While dif- 
ferent gauge choices 59ky of the Wannier states do not 
change the topological universality class of the associ- 
ated state, it can be used as variational parameters in 
the many-body ground state and can be optimized nu- 
merically by the comparison of th e Wannier ground state 
with the exact ground state. 122123 



In this section, we review the PP representation of in- 
teractions in FQH systems, and then map it to the FCI 
using the map established in ( 16 ). The PP approach was 
pioneered by Haldane in the context of FQHEj^: Suppose 
we assume all relevant degrees of freedom to be located 
in the LLL. When projected onto this subspace, which 
is a perfectly flat band, the kinetic energy is frozen out 
and the Hamiltonian thus consists only of an interaction 
term reminiscent of the Coulomb interaction 



H 



i<j 



(18) 



where the sum extends over all pairs of particles. This 
means that the interaction is fully characterized by its 
interaction energy between the particles i and j. The 
PPs are then defined by 



H = 



(19) 



where P™ projects onto a state where particles i and j 
have relative angular momentum m and is the en- 
ergy penalty for having two particles in such a state, 
taking on odd values m = 1,3,... 2Ai + 1 for fermions 
and m — 0, 2, . . . for bosons. Since the LLL is per- 
fectly flat, the m index also characterizes the distance 
between the particles. Since the interactions are of re- 
pulsive Coulomb type, we can make the V^s positive 
semidefinite by also fixing them to be real. This pseu- 
dopotential construction is particularly elegant on t he in- 
finite plane and the sphere where it was first definecP^'^, 
where translational symmetry and rotational symmetry 
are preserved and the relative angular momentum of a 
two-electron state is thus a well-defined quantum num- 
ber. In particular, this construction allows for a most ex- 
plicit connection between the Hamiltonian and the clus- 
tering properties of quantum Hall wave functions. For 
the Laughlin state wave function at = 1/3 filling, 
^ nj<j(^i~^i)^; where Zi could either represent the 
spinor coordinates on the sphere or the complex Carte- 
sian coordinates on the plane, the wave function decays 
as the 3rd power as two particles approach each other. 
This shows that two electrons are never allowed in a rela- 
tive angular momentum state of m = 1, i.e. the particles 
maximally avoid each other within a featureless liquid 
state^*. It follows that |4'l) is located in the nuUspace of 
^Laughlin = Y.t<] ^iPij- Furthermore, the Hilbert space 
at the appropriate filling only allows for one state of such 
a kind on a trivial genus manifold, implying that |^l) is 
the exact unique ground state of Hj^. Similar trial Hamil- 
tonian constructions are allowed for many other states, 
including non-Abelian quant um Hall states such as the 
Rcad-Rezayi series and otherJ^nEMI]. 

The Landau level on a cylinder can be seen as a hy- 
brid version of the spherical and planar scenario. Starting 
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from a Cartesian (a;, y) plane, we impose periodic bound- 
ary conditions along y and, in analogy to the sphere, 
quantize the pseudo-momentum along y according to the 
total magnetic flux N^, constraining the available area 
for the guiding center motion. We assume an infinite 
cylinder along the x-direction. On the cylind er g eome- 
try, the PP Hamiltonian takes the explicit fornPS 



xe 



-7 " /"^piinXi/lB piqii)i-Vj)p-iinXj/lE 



(20) 



This decomposition is valid for sufficiently short-ranged 
potentials V{ri — rj). Note that the functional form 
of Eq. [21] is different from that in Ref. |311 because we 
have used ordinary coordinates instead of guiding cen- 
ter coordinates. If one replaces all coordinates includ- 
ing their derivatives with their guiding center analogues, 
1^5'^ {ri — rj) will be replaced by the exponential tail ex- 
pression e^l''"^''^l^/^'B in Ref. EH The details of this 
calculation are shown in Appendix [K\ 

For later purposes, we also Fourier transform Eq, 
and invoke the orthogonality relation of the Laguerre 
polynomials (see also Appendix IbI) to obtain 



21 



where 7 is the aspect ratio of the cylinder. Is is the mag- 
netic length, q denotes the momentum variable along x, 
V"^ is the PP energy of a state with relative angular 
momentum m, and Lm denotes the mth Laguerre poly- 
nomial^. (Note that we have explicity included all units 
in (201 as compared to Ref. il?) Eq. 20 gives a pair in- 
teraction energy for each m. In the plane limit where 
7n — )■ kxl, the sum over n reduces to a momentum in- 
tegral along y analogous to the q integration along x, 
resulting in a two-dimensional momentum integral which 
reduces to the pair interaction 

y(r,-r,)=^y™(r,-r,) 

771 

= J2 V"'L„,{-llV^)6'{r, - r,). (21) 



'L^{kHl)V{k). 



(22) 



The above expression, which will also be rederived in 
Appendix [B] as a special case of a much more general re- 
sult obtained from first principles, enables us to compute 
the PP coefficients V^™ directly from a generic potential 
V{k). It is the starting point for the generalization to in- 
teractions involving more than two bodies, as is described 
in Section VI. 

To apply the PP decomposition to an FCI system with 
periodic boundary conditions in both directions, we have 
to compactify the open direction of the cylinder. The 
single particle states on the torus are the ^k{x, y) defined 
in Eq. 15 In this basis, the mth PP Hamiltonian has the 



matrix elements 



ni,n2,n3,n4 



= / d'rd'rj (i^„, (n)*i^„, (r,)* - ^Pn^ (n)* V'm (rj)*)[/™(n - rj){i^n, (n)i^„4 (.i-j)- i^n, (r^)^„3 (r,)). (23) 



U™' refers to a normalized PP that has nonzero projec- 
tion only in the mth relative angular momentum sector. 
They form a basis in which a generic potential V is ex- 



panded. As such, the V"^s which appear in Eq. [22] and 
other places below refer to the component of V propor- 
tional to U"\ i.e. projected onto the relative angular 
momentum sector m. For simplicity, we have denoted 
=27rn/L„ (r) as V«(r), n = 1,2, L^Ly. 
We can use the map (16) defined in Sec. Il]to define 
the corresponding PP Hamiltonian in FCI, which has the 
second-quantized form 

»l>n2,n3."4 

n+;iez,n+'2ez 



Here, 



(25) 



^. The matrix element U\ 



fied to the form of UV 



is simpli- 

— rrm j _ i„ 

~ '^n+li,n-li,n~l2,n+l2 

the magnetic translation symmetry. (More discussion on 
the magnetic translation symmetry will be presented in 
Sec. |V]) Depending on whether we consider the torus 
or cylinder geometry, the sites along the main cylin- 
der axis labelled by n are assumed to obey periodic or 
open boundary condition, respectively. For the cylin- 
drical case, Eq. [24] can be brought into a bosonic pair 
creation form given by C/;™^ = (7k^6™6™, where 



p=0 



(j-p)!(2p+l)! 



so that 



(27) 



where 5™ = E„+/ez ^ra"-/an+i- Here, g ^ JY»,2 is a 



is the annihilation operator of the Wannier state constant with units of energy and k 



(27r)3/2 
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has been set to ^ lattice constants in the latter equality 
in accordance to Ref. 1211 In the following, Is will be 
expressed in units of the lattice constant unless it appears 
in the combination /^V^ or l%q^, where g is a momentum 



variable. The complete derivation of Eq. 27 can be found 
in Appendix [C| 



The Hamiltonian in (24) will be the starting point of 



Sec. |IV| when we expand different FCI models into this 
PP form. Its m = 1 case has been previously used to de- 
fine low-dimensional Mott-type models with bare onsite 
hardcore potentials at fractional fiUing"^^ ^31. The PP 
on the torus can be found by summing over all the peri- 
odic images of V^™^ (referred to as in Appendix |d|) 



satisfying li + I2 mod 2LxLy — 0. This constraint can 
be generalized to the case with more than two bodies, as 
shown in Appendix [E} 

For finite-size investigations on the cylinder or the 
torus, we have to keep in mind that relative angular mo- 
mentum is no longer a well-defined quantum number, as 
opposed to the case of the sphere or the plane. The 
parameter m in (21 ), which corresponds to the exact rel- 



ative angular momentum as we take the planar limit, 
determines the order of the derivative acting on the hard- 
core potential via the degree of the Laguerre polynomial. 
This corresponds to a Taylor expansion of the interaction 
in momentum spac^^. While its interpretation as the ex- 
act relative angular momentum is absent, it can still be 
employed as an expansion parameter for short-range in- 
teractions on a sufficiently large torus or cylinder. To see 
this in terms of the Hilbert space basis, we describe rela- 
tive motion states on the torus by relative motion states 
on the plane. The latter can be exactly classified via the 
relative angular momentum m which is proportional to 
the interparticle distance in that relative state r^. For 
fm/Lx,rm/Ly « 1, the overlap of the torus and planar 
relative motion states goes to unity, effectively reestab- 
lishing the notion of torus relative angular momentum 
for short distances. Still, this approximation becomes 
invalid for higher values of relative angular momentum. 
At the Hamiltonian level, this is reflected by the over- 
completeness of the PPs V^. This occurs because the 
interparticle distance ^ Lx,Ly that characterizes a 
relative angular momentum state will no longer be well- 
defined when rm is comparable to the system lengths of 
the torus. A quantitative treatment of the overcomplete- 
ness bounds can be found in Appendix [P] 

A deep insight to note is that even though an ex- 
act angular momentum quantum number cannot be de- 
fined, the clustering property of the quantum Hall-type 
wave functions are still fixed appropriately at these finite 
size manifolds such that they can be exactly located in 
nullspaces of PP Hamiltonians. This was elegantly shown 
for the torus by Haldane and Rezayi^^, which we illus- 
trate for the Vi PP at v — 1/3 filling: demanding that 
the many-particle ground state pays no energy due to 
Vi , it necessitates that the wave function decays to third 
power as the particles approach each other. Oddness due 
to fermionic statistics and boundary conditions on the 



torus automatically restricts the functional form of the 
wavefunction to be ^E* ~ Y\i<j'^i^i ~ ^ik)'^: where d is 
the odd Jacobi theta function. The groundstate is thus 
forced to be of Laughlin type. This fixes N^p — 3 zeroes 
of the wave function, where the remainder 3 constitute 
the topological center of mass degeneracy of the Laugh- 
lin stat^35]_ A similar discussion can be pursued for the 
cylinder, where the center of mass degeneracy is absent 
but the clustering property of the wave function leads to 
the same finding for the remainder functional form of the 
wave functional. 

All in all, the properties of the pseudopotential expan- 
sion sets the stage for numerical investigations of short- 
ranged interactions of FCIs as well as the defining of 
trial Hamiltonians for new quantum Hall-type fractional 
Chern phases, both of which will be pursued in the fol- 
lowing. 



IV. PSEUDOPOTENTIAL EXPANSION OF 
FRACTIONAL CHERN INSULATORS 

A. Model Hamiltonians 

We apply the PP expansion to two model FCI Hamilto- 
nians, the checkerboard (CB) model introduced in Ref.lTTl 
and the honeycomb (HC) model introduced in Refs.lHland 
l47l Both models possess an almost flat (dispersionless) 
band which mimics the LLL in a FQH system. In an 
FQH system, the Coulombic electron interactions lift the 
macroscopic degeneracy of the LLL, leading to a topolog- 
ically degenerate groundstate. In the same spirit, we add 
Hubbard-type interaction terms to our model Hamiltoni- 
ans such that 

H = Hq + Hint , 
Hint = ^Hnn + (1 — X)Hnnn 



— Hq + ho \ X ^ niUj + (1 - A) 



(('y)> 



(28) 



where A characterizes the relative strengths of the nearest 
neighbor (NN) and next-nearest-neighbor (NNN) inter- 
action terms, h^, a parameter with units of energy, sets 
the overall magnitude of Hint- The single-particle term 
Ho gives rise to the almost flat band and provides infor- 
mation for the construction of the Wannier basis. This is 
the basis we will use for expressing the two-body interact- 
ing term _ffi,it in the same basis as the PP Hamiltonians 
U^,U^, etc. as denoted in (p4). 



We consider interactions that are much larger than the 
bandwidth of the almost flat band, but much smaller 
than the interband gap. In this limit, for a partially filled 
flat band, we can ignore the coupling to the upper band 
and only study the effect of interactions in the subspace 
of the flat band. With this picture in mind, we expand 
Hint in terms of the PPs in the Wannier basis of the 
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partially filled band. As long as the bandwidth of the 
almost flat band is much smaller than the interaction 
strength, we can ignore the bandwidth and consider only 
the interaction term. 

The checkerboard (CB) lattice model consists of two 
interlocking square lattices displaced (1/2,-1/2) sites 
relative to each other (Fig. [S]) . Its noninteracting Hamil- 
tonian H^^ consists of NN, NNN and NNNN hopping in- 
teractions parametrized by hopping strengths t, t' and t" 
respectiveljB^. The NN interaction exists between sites 
belonging to different sublattices and carries a phase 
giving rise to the time-reversal symmetry breaking nec- 
essary for a nonzero Chern number. Both the NN and 
NNNN interactions exist between different sublattices, 
leading to off-diagonal terms in the noninteracting Hamil- 
tonian. Explicitly, in sublattice space. 



FIG. 3. (Color online) Lattice structure of the CB(left) and 
HC (right) models. Sites colored differently belong to differ- 
ent sublattices. The unit cells are demarcated in red. For the 
CB model, NN interactions are between different sublattices 
while NNN and NNNN interactions occur within the same 
sublattice. For the HC model, both the NN and NNNN in- 
teractions occur between different sublattices, but the NNN 
interactions act within the same sublattice. 



where 



k k 
di = —4:t cos (p cos cos , 



0,2 = — 4t sm (j) sm — sm , 



in the same sublattice. NNNN interactions which occur 
for diametral sites on the same hexagon involve different 
sublattices (Fig. [s]). After an analogous gauge transfor- 
mation. 



where 



di = —t(l+cos kx+cos ky)—t'' {cofi{kx+ky)+2 cos{ka:—ky)), 



dj — —2t' {coskx ~ cosfcy). 

The expression for dp is irrelevant because it is not needed 
for the computation of the Wannier basis. We set t = 1, 
i! = 1/(2 + y/2) and = 7r/4 as in Ref. [T7|to maximize 
the flatness of the bottom band. We can explicitly see 
why a nonzero 4> is necessary for having a topologically 
nontrivial gapped model: as the Chern number is given 
by Ci = ^ / d^ATd • {d^d x dyd), it can only be nonzero 
if none of the di's is identically zero. 

Notice that Hq^ is not of Bloch form since the d^'s 
do not obey the periodicity of 2n. This is because some 
sites are noninteger lattice spacings away from each other 
(Fig. [3]). We can remedy this by shifting one sublat- 
tice site on top of the other within a unit cell. Mathe- 
matically, this corresponds to a gauge transformation of 
"^Ib ~^ c|,^e~'(*''="^*''«)/^ where B refers to one of the sites 
within the sublattice. After the gauge transformation, 

di — —2t[cos(l)+cos{kr^+ky—(j))+cos{kr^+(t))+cos{ky+(f))], 
d2 = —2t[— sin (p+sm.{kx+ky—(l>)+sm.{kx+(f>)+sm(ky+(j>)], 

^3 — ~2t'{coskx — cos ky). 

The noninteracting part of the honeycomb model is de- 
fined similarly. The unit cell consists of two adjacent 
sites. The phase cj) is carried between NNN sites, which lie 



^2 = i(sin kx + sin ky) + t" sin(fc2, -I- ky)^ 

d^ — 2t'sin(/)(cosfc2: -l-cosfcj, -t- cos(fca; — ky)). 

The values for the NN, NNN, and NNNN hoppings are 
given by t = 1, t' = 0.6, and t" = -0.58, = 0.47r such 
that the flatness of the band is optimized'*^. We stress 
that the optimization of these flatband parameters is not 
necessary for performing the PP expansion, while it is 
physically relevant to increase the stability of any FQH 
state present in the system. 

B. Expressing Hint in the Wannier basis 

We have the necessary ingredients for expressing the 
interaction part of our model Hamiltonians in the same 
basis as the PPs in FQHE. First, we can perform a 
Fourier transform on such that it becomes 

Hint = A/io ^ rLiirij + (1 - A)/io ^ niUj 
(■y> ((ij)) 
= ^T/°^(g)n_,,,n,,, (29) 

where q is an internal momentum variable, a, h are the 
sublattice indices, and n^a = Efc.fe, 4+9,aCfca- V''^{q) 
denotes the gth Fourier component of the interaction be- 
tween the sublattice index a and h. This is an expression 
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quartic in the creation and annihilation operators c, in where the normahzation factors Nqa-^^a^x^x^k follow 

the momentum/sublattice basis. Since we are only con- from (30). We obtain 

sidering interactions within the flat band, we project out 

the upper band and only keep the matrix elements of 

Uqa in the flat band. After projection, the annihilation 

operator Cka can be expanded in the Wannier state basis: 



ka^^y-"- W^27rn/L„) a« + upper band. (30) 



In this representation, the density operator becomes 

nqa = ^ I C^+^^a^fca 



Er-1/2 TT* t 

X / , '^nixi.k^+q^,a,ky+qy"-xi,ky+qy,ni 
k^ky n\xi 

^ ^ Un2X2k^aky^X2kyn2 

n2X2 

^ 5Z '^<iaia2XiX2kyO'l;-^^ky+qy,ai'^X2kya2^ 

(31) 



kyaia2XiX2 



Hint = KT'^-g.a? 



^ ^ ^ ^ ^qy — k2y-\-kiySk2y — kiy—ksy + k4y^—qxiaia2X2k2y-^qa3a4XsX4k4y(^\c^0,^ 

qx {xj},{Kjy},aj 

KiK2K3Ki 

- X! X! 'l' ^2)(a(„-ii)/2a(n+/i)/2)^(a(n'-;2)/2«(n'+Z2)/2), 

n'n I1I2 



(32) 



where Ki = kiy + 27rxiCi 



kiy + 2ttxi. (The other 
KiS are defined analogously.) A quadratic term has been 
dropped in the final step because it can be absorbed 
into the noninteracting part of the Hamiltonian. Due 
to fermionic statistics of the c„±; operators, we can an- 
tisymmetrize iJint, leading to the matrix elements 



h{n,n',lij2) = f{n,n',li,l2) - /(n, n', ^2) 

-f{n, n', h.-h) + /(n, n', -h, -h), 



(33) 



which are manifestly antisymmetric in li and I2, just like 
the PPs U{^i^. We see from Eq. [32] that h{n,n' ^h.h) 
corresponds to a pair hopping interaction on a line. Two 
particles with the CM "position" n' separated by I2 sites 
simultaneously hop onto new positions with CM position 
n separated by li sites (see also Fig. [6]). More discussions 
on the physical interpretation of this interaction will be 
presented in Section [V} 



C. Pseudopotential expansion of the interaction 
term 

While the PP matrix elements UJ^i^ depend only on 
li and ^2, the FCI interaction Hamiltonian matrix ele- 
ments in the Wannier basis /i(ni, 712, /i, ^2) also depend 
on Til and 712. Hence only a part of h{ni, rt2, h) can be 
expanded in terms of PPs. This important fact can be 
understood in terms of magnetic translation (MT) sym- 
metry breaking, which will be analysed in depth in the 
next section. Here, we shall concern ourselves with the 
terms that can be expanded in PPs, defined by 



Hp{n,n\li,l2) = 



2Ly 



^ y N=i 



h{N,NMM)- (34) 



Hp vanishes for n ^ n' and does not depend on n when 
ni — n2 — n, as required. The sum runs from 1 to 2Ly 
because h{n, n, Zi, ^2) is periodic in n with period 2L^Lj,, 
as evident from the periodicity of a(n±i)/2 in Eq. 32 
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We would like to expand Hp in an orthonormal basis 
of PPs C/^, , etc. However, this expansion is only 
unique and thus meaningful if we include PPs with m 
bounded by a certain Wmax- This is because the inclusion 
of higher PPs can yield an overcomplete operator basis, 
a consequence of the finite size of the torus geometry 
explained in Section |III[ The truncated PP basis is no 
longer complete, but we can still perform a PP expansion 
of Hp by writing 



Hr. 



H^ =: H 



pscudo 



(35) 



and finding the coefficients aj that maximizes the nor- 
malized overlap {Hp, H pscudo) ■ The overlap is taken by 
summing over all 1^21 < LxLy since Hp has a period 
of 2LxLy. Specifically, for any two Hamiltonians H and 
H' that respect MT symmetry, 



{H,H') = 



Hi.uH' 



hh 



^hh ^?il2 ^hh ^hl 



(36) 



The term i?> consists of the part of i?int that does not 
break MT symmetry, but still cannot be uniquely ex- 
pressed in terms of the PPs. It includes, for instance, 
hoppings that occur over lengths comparable to the size 
of the torus. 

When the [/™s form an orthonormal basis, the 
coefficients a, that maximize the normalized overlap 



[Hp, iJ] 



pSCUdOy 



can be determined as 



{Hp,U"^). 



From Fig. |4] we see that the percentage of Hp that can 



be expanded as PPs X]j=o°^ I^ T has the maximum value 
for 0.94 for the CB NNN interaction (Fig. As ex- 
pected, the first PP has the highest weight, which favors 
the possibility of simple FQH states such as Laughlin 
states. With the relative angular momentum m being 
proportional to interparticle distance, Um is expected to 
decay faster with |/i|,|Z2| as m increases. This will be 
shown in more detail in Appendix [C] It is notable, how- 
ever, that the second neighbor coupling leads to a bet- 
ter overlap with the first PP than the nearest neighbor 
Hamiltonian. This indicates that the PP Hamiltonians 
mapped to FCI systems are not simple density-density 
interactions and their matrix elements in real space lat- 
tice site basis can exhibit a nonmonotonic dependence on 
distance. More specifically, this is because U^^i^ ~ ^1/2 is 
not strongly peaked around li — ±^2 in h — I2 space, as 
shown in Fig. [5] unlike the NN interaction. Since li — ±^2 
corresponds to g = (as defined in Eq. 32), we see that 
the NN terms are "too local" for a good overlap with U^. 



In general, the matrix elements U'l 



hh 



{hh)"", so U" 





FIG. 4. (Color online) The pseudopotential expansion of the 
checkerboard (upper) and the honeycomb (lower) models for 
= Ly = 6. The normalized overlap aj is plotted against 
A, where -ffint ~ \Hnn + (1 — X)Hnnn- We see that the NN 
and NNN terms of the CB model exhibit marked diferences 
in their pseudopotential expansions, with the MT symmetry 
conserving part of its NNN term consisting almost exclusively 
of the and terms. This can be understood by studying 
its distribution of matrix elements (Fig. [5|. 



V. THE EFFECT OF MAGNETIC 
TRANSLATION SYMMETRY BREAKING 

In this section, we analyze the origin of the terms in the 
FCI Hamiltonian that cannot be expanded into PPs. We 
review how magnetic translation (MT) symmetry in FQH 
system constrains the form of its two-body interaction 
terms, and investigate how this picture is generalized to 
the FCI case. 



A. Origin of Magnetic Translation Symmetry 
breaking 

Consider an L^Ib x Lyls torus geometry which has 
been discussed in Sec. II. In the Landau gauge = 



becomes more localized at li = ±^2 for higher m. 



0, A, 



Bx, the covariant momentum operators are 



11 




FIG. 5. (Color online) Density plots of the NN (left) and 
NNN (center) terms of \Hp\ for the CB model. The horizontal 
and vertical axes represent h — I2 and h + h respectively. 
Regions with relatively large \Hp\ are colored red. (right) 
is plotted for comparison. The NNN term evidently bears 
more resemblance to . 



Px = —idx, Py = —idy—Ay which satisifes [Px, Py] — iB. 
The Hamiltonian -ff = 5™ i^x + Py) has two translation 



symmetries , Ty defined by 



o— • 



-o-»-«-o— o-»-«-o- 

V> VJ^ VJ^ 



FIG. 6. (Color online) The two types of MT breaking hop- 
pings. Top: A hopping process that changes the CM position. 
Bottom: Two hopping processes that preserve the CM posi- 
tion but break magnetic translation symmetry because the 
hopping at different CM position n (solid and hollow arrows) 
have different amplitude. 



1^ = e e 



TB - JPyTi 

±y — C 



(37) 



Ty is an ordinary translation while is a transla- 
tion in X direction by Is/Ly accompanied by a gauge 
transformation. The translation can only be defined in 
units of Is/Ly so that the change of gauge potential 
Ay = Bx B{x + Is/Ly) can be cancelled by a gauge 
transformation. The action of , Ty on the basis wave- 



functions (15) is 



Kl 



JK/L, 



(38) 



For a general two-body interaction with iJjnt in the 
form of — V r/nin2n3«4„ „ 

the condition [T^,Hi^t] = requires ni + n2 = 

«3 + «4, since (T^f) ^ a^^a^^an^amT^ = 

4i42a«3a»4e^'''^"'^"'""'~"''^^'^"'^''- The con- 

dition [Tf,iJint] = requires t/»i»2"3n4 = 
^ni+i,«2+i,n3+i,«4+i Therefore, the magnetic trans- 
lation symmetry T^f and Ty determines the CM 
conservation (ni + n2 — nj, -t- 71.4 or n = n') and (one- 
dimensional) translation symmetry of the interaction 
Hamiltonian in FQH states, i.e., n-independence of the 
interaction matrix elements. 

By comparison, in the lattice model, we only have the 
lattice translation symmeties which commute with each 
other. The action of the lattice translation T^, Ty acts on 
the Wannier basis as 



Ty \Wk) = e^^ \Wk) 
Tx \Wk) = \Wk+2.) ■ 



(39) 



the matrix element of two-body interaction ^"i"2«3n4 
satisfy 

Jjnin2n3n4, _ jjni+Ly ,n2 + Ly ,n3+Ly ,n4+Ly 

j^«in2n3"4 ^ if ni -t- 712 7^ "3 + ^4 mod Ly. (40) 

The magnetic translation symmetry breaking in the lat- 
tice models (Fig. [6]) is also related to the non-uniform 
Berry curvature in momentum space. As was discussed 
in Sec. |llj the CM position of the Wannier state \Wk) is 
determined by the flux of the Berry's phase gauge field 
Px{ky). If the system has magnetic translation symme- 
try, \Wk) and \WK+2'n/Ly) are related by T^ , so that 
Px{ky) must depend on ky linearly. As a result, we ex- 
pect MT symmetry breaking whenever the Berry curva- 
ture is nonuniform in momentum space, which is the case 
in a generic CI. 

In addition, MT symmetry breaking will still be 
present even in the hypothetical case of perfectly flat 
Berry curvature. This is because the Wannier basis is 



not perfectly local. Recall from ( 32 ) that 



ni + n2 = n — kiy + + (a^i + X2)L 



^y, 



713 -I- 714 = 7l' = kly + k2y + {x-i + X4)Ly 



(41) 



(42) 



where the XiS are the lattice sites of the original i/int- 
CM nonconserving terms occur where X1+X2 7^ 2^3 +2:4, 
when n and 7?,' differ by a multiple of Ly. These terms do 
not appear in the original real-space basis where i?int oc 
'^ij^i'Tij — —^ij c^iC^jCiCj annihilates and creates two 
particles at the same position. However, our Wannier 
basis functions generically have exponentially decaying 
tails on both sides of their peak x, which produce CM 
nonconserving and thus MT breaking contributions. 



Comparing Eq. |39] with Eq. |38[ we see that in the map- 
ping from FCI to FQH defined in Sec.|Ilj Tx,Ty is mapped 

to (Tf )^'' and {T^)^\ respectively. Therefore, in the 
lattice model, the translation symmetries only require 



B. Numerical results on MT symmetry breaking 

We present the numerical results on MT symmetry 
breaking in our model Hamiltonians. Define the resid- 
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n 

FIG. 7. (Color online) A plot of Sl for = Ly ^ 6 for dif- 
ferent model Hamiltonians. As evident from the dominance 
of the peak at n = 0, the amount of MT symmetry breaking 
largely stems from CM conserving terms. There is little dif- 
ference between the degree of MT symmetry breaking in the 
different models. 



ual 

Hycs{n,n',li,l2) = h{n,n',li,l2) - ^p(^i,^2) 

where, as before, h{n,n' ,li,l2) denotes the FCI interac- 
tion Hamiltonian expressed in the Wannier basis. i?ics is 
the part of the FCI interaction Hamiltonian h{n, n' , li, I2) 
which does not satisfy MT symmetry required by the 
PPs. Obviously, H^^s = if /i is one of the PPs, since h 
will then be equal to Hp. The quantity 

allows us to track the origin of MT nonconservation. (5q 
comprises the elements of H-^cs satisfying n — n' . As 
defined in Eq. [34j these are the elements which are in- 
dependent of n. (5q hence represents the fraction of ma- 
trix elements that are CM conserving but MT symmetry 
breaking. For n 7^ 0, 5^ represents MT nonconserving 
contributions that likewise do not respect CM conserva- 
tion, is plotted in Fig. [t] for various model Hamiltoni- 
ans, for a system size = Ly = 6. The results remain 
almost unchanged when L^ and Ly are varied as long as 
Lx ~ Ly > 3. 

From the enhanced peak at n = 0, we conclude that 
most of the MT symmetry breaking occurs when CM 
is conserved. This happens because our maximally lo- 
calized Wannier functions (WFs) are still mostly peaked 
at one site. The subdominant contributions from (5^^ for 
n = ±1 can be attributed to the finite tails of the WFs 
one site away from their center of mass. Indeed, S'^ be- 
comes exponentially small for |n| > 1. While the over- 
all extent of MT symmetry breaking originates from the 
nonuniformity of the Berry curvature, its relative contri- 
bution to Sf^ for different n is dictated by the localization 
properties of the WFs. 



C. Discussion 

The decomposition of the FCI Hamiltonian into pseu- 
dopotentials is only exact in the thermodynamic LLL 
limit of zero bandwidth and homogeneous Berry curva- 
ture. For the generic model, the FCI Hamiltonian can 
only be partly decomposed into pseudopotentials, which 
we then discuss along general FQHE pseudopotentials 
on the cylinder or torus. From our calculations, the de- 
viations are significant, suggesting that at least for the 
spectrum above the elementary low energy quasiparticle 
regime, there is no clear similarity between FCI and FQH 
systems. However, entanglement signatures of incom- 
pressible liquid phases, such as the entanglement spec- 
truirP^ with the emergence of an entanglement gajP^, 
show strong similarities of FCI ground states to their 
FQH analogue s, eve n in terms of the counting rule of 
low-lying state42U33)_ This is astonishing from the view- 
point of PPs, as the FCI and FQH models at the bare 
level could only possibly agree to the extent of the PP 
decomposable components of the FCI Hamiltonians. 

Such an apparent discrepancy between the points of 
view of Hamiltonian and entanglement measures can be 
interpreted as a consequence of a renormalization group 
flow. As high energy modes are integrated out, the low 
energy physics of FCIs supposedly flows towards FQHE 
type scenarios, with reemergent symmetries such as mag- 
netic translation group which is conserved in the FQHE 
but broken in the FCI at a bare level. Recast into pseu- 
dopotentials, it suggests that the PP non-decomposable 
part of the FCI Hamiltonian at the bare level should 
decrease upon renormalization, while the ratios of pseu- 
dopotentials of the PP-decomposable part of the bare 
interactions might deviate from the PP ratios in the low 
energy theory. This implies that even if two FCI Hamil- 
tonians have similar PP ratios at the bare level, they can 
still differ considerably in their low energy description, 
and hence their propensity to host FQHE-type incom- 
pressible states (Fig.[7|. This interpretation is consistent 
with the common theme from FQHE numerical studies 
that the data quality of entanglement spectra and its 
characterization of the bulk and edge mode properties is 
not significantly correlated to the spectral sharpness of 
the Hamiltonian spectrum, and partly anticipates the en- 
ergy spectral flow^. Ultimately, only the joint confirma- 
tion of both entanglement and Hamiltonian measures will 
justify true evidence for a fractional topologically ordered 
phase in the FCI models. From a low energy perspective, 
the seemingly clean finding from entanglement measures 
might not yet rule out that the inhomogeneous berry cur- 
vature induces a fiow to a liquid different from FQHE, as 
may also be seen by hints such that the hierarchy liquid 
construction cannot be established for the FCIs as in the 
FQHE cas^il From the perspective of energetics, the 
effective pseudopotential weights of the FCI models in a 
low-energy theory are likely to be strongly modified due 
to "integrating out" the PP non-decomposable part of 
the bare model, which can also provide an explanation 
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for the parameter trends of the stabihty of FCI phases as 
a function of system parameters'^ . 

VI. MANY-BODY PSEUDOPOTENTIALS 

A. Construction of many-body PP trial 
Hamiltonians 

While we have only explored two-body interactions so 
far, PP expansions are also well-suited for many-body 
interactions. From the established knowledge in FQHE 
systems, it follows that various interesting FQAH liquids 
are located in the nuUspace of certain many-body PPs. 
In theory, we can construct many-body FCI Hamiltoni- 
ans that exhi bit Pfafhan , Read-Rezayi etc. groundstates 
from such pppMEHlSl. 

The first task is to generalize Haldane's PPs for FCI 
models to more than two-body interaction^^. For two 
particles in the LLL and a translationally invariant po- 
tential V{k) = / e~'^'^'^V{r)dr, the mth Haldane pseu- 
dopotential coefficient is given by Eq. [22] 

V- = inll J ^^e-^l'^\^{lie)V{k) (43) 

so that the potential with c>c 5mn is given by U"''{k) = 
VoL„i{lgk^). As before, Vq is a constant with units of 
energy. In the plane limit, the U"^s form an orthogonal 
basis which one can use to expand a generic short-range 
potential. 

When an interaction involves more than two parti- 
cles, additional complications arise. To begin with, there 
are different ways of choosing to assign relative distance 
variables, or, angular momentum. (For two-body inter- 
actions, there is a unique assignment, as one degree of 
freedom drops out due the CM conservation.) When 
there are 3 particles, only one degree of freedom is elim- 
inated due to CM conservation. As such, an ambiguity 
remains in choosing the many-body analog of relative 
angular momentum. This ambiguity is mathematically 
manifest when one tries to generalize Eq. [43] In the case 
of 3-body interactions, there will be two momenta ki 
and /c2 integrals in the above expression, and one has 
to chose the new expression to involve Lm(fci), ^771(^2)' 
Lm{{ki — fe)^), or a combination of these. 



This formal ambiguity similar to coupling multiple an- 
gular momenta does not induce physical complications 
as we formulate generalized Haldane pseudopotentials 
(GHPs). Its detailed first-principle derivation can be 
found in Appendix [B] We constrain ourselves to the ap- 
plication of GHPs to the total relative angular momen- 
tum iV-body PP Hamiltonians: 

c-w^i„(5|i^)^ m 

Here, U™'{k) is the N-body interaction potential that has 
a total relative angular momentum of m, with k being the 
momentum conjugate to the total relative coordinate w 
defined by 

- - ]^ E - Zn) = ZN. (45) 

n=l 

where Zi = Xi — iyi are the complex coordinates of the N 
particles. 

In real space, the pseudopotential U"''{w) cx 
J dke'^^''^U"^{k) depends explicitly on the positions of 
each of the N particles. If we select one of the N parti- 
cles, the total relative angular momentum is the sum of 
the relative angular momenta of the other iV — 1 particles 
relative to the ffist one. Indeed, we see from Eq. [45[that 
w represents the relative seperation between particle N 
and the CM of the rest of the particles. Note the appear- 
ance of the factor 2(n-i) ' which is essential in obtaining 
the correct expressions for the PPs. (It will be derived 
in detail in Appendix [b] ) 

The next step is to express Um{k) in the LLL Lan- 
dau gauge basis. Since the latter is the analogue of the 
Wannier basis, once we have done so, we are able to read 
off the Wannier basis matrix elements of a many-body 
FCI Hamiltonian exhibiting Pfaffian or other more ex- 
otic groundstates. Here we shall perform this explicitly 
for iV = 3 bosons. The computation with fermions or 
more bodies is conceptually similar. According to Eq.[44] 
we rescale the magnetic length P by |P. Hence, U™ be- 
comes yo-^m(|fc^), where k is the momentum conjugate 
to the total relative coordinate ((?^i — r2) -I- (ra — r2))/2. 
In the basis of LLL Landau gauge eigenfunctions. 



jjm 

«'^''l'^^^2d^r3V^1,^(ri)V'i,(r2)V'Jj3(?-3)im (^|^|Vf(^^+^3)/2_^^]^ V'„4(^3)V'„5(r2)V'„6(^l) 

«E / lf[d'nL„,(J{p~qri%^ e^'-^'-^~'-^^e^^-('--'-3)<(ri)<(r2)Vi3(^3)^„,(r3)V'„,(^2)^„Jri). 

(46) 
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The sum refers to a symmetric (antisymmetric) sum 
over all permutations a assuming the particles are bosons 
(fermions). We have k = q — p because 



gig-(ri-r2)gjp-(r2-r3) 



As before, ipnir) 



gi(g-p)-((ri+r3)/2-r2)gi(p+<3)-(ri-r3)/2 



(47) 



tL/f 



and K, — 



is a dimensionless ratio that is small in the limit of 
large magnetic fields. 



convenient form. As its computation is instructive for the 
generalization to more complicated cases, we explicate it 
in Appendix [e] The 3-body bosonic PP C/" which hosts 
the Pfaffian ground state is given by 



where 



E 

ni=3RmodN 



E 



^ni *-"n-2 



E 

ni +712+^3 —3i?mod A'' 



E' 

s,t 



'-Wst 



L-ni '-'712 ) 



(48) 



in Ref. f5S! However, the GHP formalism here can also 
generate higher PPs of multi-body interactions, and this 
result is just its simplest case. Note that the summation 
constraint ni + n2 + = 3i?mod7V can also be imple- 
mented as constraints over s,t, as in the 2-body case 
described in Ref. |43l 

We see that is positive semidefinite. Since the Pfaf- 
fian state |Pf) resides in its kernei^^^, we have 



The integrals related to (CI) can be simplified to a ^s well as 



(Pf|C/°|Pf) = 



fe^lPf) 0. 



As a consequence, the Pfaffian state is not just the ground 
state of C/° but is also annihilated by all b^s. Analo- 
gous to the two-body case, U° consists of all three-body 
processes that conserve the CM R. The physical posi- 
tions of the particles relative to their CM is described 
by i? — {rii + Nsi), where Ui is the real space basis in- 
dex before enforcing the periodicity. When k is suffi- 
ciently small, contributions from periodic images can be 
discarded, and we are left with a simpler expression with 
\R-{n, + Ns,)\ <N = L^Ly. 
Using our setup, we can likewise calculate higher PPs. 
vanishes due to the symmetric exchange statistics of 
the bosons. From Appendix |E) we have 



with Wst = Ej"-'-* - T,i<jK'riji "i = "-i + sN, n'^ = 
n2 + tN, and n(j = na - N{s + t). The N = Lr,Ly pe- 
riod originates from the periodicity of the Wannier basis 
on the torus. This result has been previously obtained 



Y.'b^RbR. 



where 



bR 



E 



El 



2k' 



W.^t e r 



■-•ni "--712 ; 



(49) 



with Wst = "Ylii '^'i ^ "^iKj "-'i = + s-^i "2 = ^2+ tN, and ng = 713 — N{s + t) as before. 

I 



B. Discussion 



Generically, the FCI models will only contain two-body 
interactions, such that our many-body pseudopotentials 
are not for primary use as an expansion tool for bare 
Hamiltonians. This situation is similar for the FQHE 
case. There, the many-body PPs are primarily used to 
provide trial Hamiltonians for which various FQHE states 
are exact null modes. As stated before in Section III, the 
exact null mode property is tied to the clustering prop- 
erty which is enfored by the A'^-body trial Hamiltonian 
as N particles approach each other. This will be an in- 
teresting analogous application for the FCI scenario. By 
use of our many-body PPs, various trial states can be re- 



alized in a FCI model, including fractional Abelian states 
such as the Laughlin series, fractional non- Abelian states 
such as the Read-Rezayi series, but also even more exotic 
states such as a FCI Gaffnian statePJI^. For the Pfaf- 
fian FCI state, for example, it is interesting to study the 
interpolation of the exact FCI Hamiltonian to a generic 
FCI model and see whether adiabatic connectivity can be 
reached both at the level of Hamiltonians and entangle- 
ment spectr£p2i. In addition, as soon as interband scat- 
tering would be considered, effective many-body interac- 
tions are naturally generated within generic FCI models, 
which provides yet another area where many-body FCI 
PPs could be employed. 

Furthermore, following the construction introduced in 
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Ref. |3H1 for the hardcore potential, we can map all these 
states to effectively one-dimensional models (see also Ap- 
pendix|B]and[E]), and even higher dimensional generaliza- 
tions-^ , of featureless Mott insulators with exotic ground 
state and quasiparticle properties, which deserve investi- 
gation in their own right. 



VII. CONCLUSION 

We have developed a pseudopotential formalism for 
fractional Chern insulators. Starting from the FCI Wan- 
nier state representation, we have employed the FQHE 
formalism on the cylinder and torus to define two-body as 
well as many-body PPs. We have decomposed bare FCI 
models into PPs and find that generic situations give rise 
to PP-decomposable and PP-nondecomposable parts of 
the FCI Hamiltonian, and discussed their interplay with 
partial breaking of magnetic translation group, which ap- 
pears to be reemergent at low energies in the FCI models. 
We have defined many-body PPs to establish a basis for 
studying further unconventional FCI phases via appro- 
priately designed trial Hamiltonians. We believe that 
it will be interesting to employ the PP perspective to 
provide a complementary tool to entanglement measures 
and to further develop our understanding of fractional 
Chern phases, as well as to extend i ts applic ability to 
Chern bands of higher Chern numb eiEsMHMl and two- 
dimensional fractional topological insulator model^^^^^^. 
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Appendix A: Equivalence of the Landau level and 
Cartesian description on the cylinder 

We reconcile the two expressions for the delta function 
potentials on the cylinder appearing in Refs. |38J and IHl 
The following arguments are valid for generic Ib, which 
has been set to unity in the following. We show that 

{0\6^ir - r')|0) = p-('-''-'-p)'/2^ (Al) 

where |0) ^ /(z)e~l^l^/^ denotes the LLL wavefunction 
and Tp = {xp,yp) are the guiding-center coordinates de- 
fined by 

Xp • + 7^ -7 

Vp -■ -lOz + Y : 

with z = x — iy. These are the LLL-projected coordinates 
since Zp = z, Zp =: 2dz Here, the normal-ordering 
symbols ":" indicate that any derivative within them 
does not act on the Gaussian factor e"'^' present in LL 
wavefunctions. We also note in passing that [Xp,yp] — i, 
i.e. that the guiding-center coordinates do not commute. 
Consider a general interaction 

Let us project it onto the LLL |0) = f{z)e~^^^^/^, i.e. 

{0\Vir-r')\0)^ J ^nfc)(0|e*('-''')|0). 

The quantity in the angle brackets can be evaluated in 
an analogous way as Eq. |B10[ Specifically, we have 

= e-l^-|'/2(0|e7l(^^'+'='')|0), (A2) 

where k = kx — iky and a, b, , are lowering and rais- 
ing operators of angular momentum and LL (see also 
Appendix [b|) . Due to the specific form of |0), we also 
have 



6t|0> = ^(^ - 29,)/(z)eH-lV4 = : z : fiz)e-\^\'^' ^ -^Zp\0). (A3) 



Hence, If we start from a delta function in real space, V{k) — 1, 

(Ole^'^'nO) = e-l'=l'/2(o|et('=^='+'=^)|0) 

= e-l'=l'/2(0|e^'="-''|0). (A4) 
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we obtain 



Allowed coordinate transforms Ri 



(0|^2(^_/)|0) 



g-|fe|V2gjfe-rp 



V2 



Higher PPs in the LLL will thus be of the form 
Lm(Vp)e^''''''~''3'-' To summarize, the LLL projection 
leads to two types of modifications. Firstly, the delta 
function is replaced by a Gaussian in guiding-center co- 
ordinates. Secondly, the derivative will also be taken with 
respect to guiding-center coordinates. 



Appendix B: First-principle derivation of 
Generalized Haldane Pseudopotentials 

We generalize the Haldane pseudopotentials to iV > 2 
bodied interactions, i.e. GHPs, in fractional Chern 
bands. The main calculational steps in this appendix in- 
cludes (i) a derivation of the GHP starting from an orig- 
inal interaction Hamiltonian involving arbitrarily many 
bodies (ii) a clarification on how the total relative an- 
gular momentum can be defined through an appropriate 
change of coordinates. The results of this appendix will 
be directly utilized in Appendix |E] for explicit calcula- 
tions of the first 3-body PPs , which always take the 
factorized form ^ X^n ^n^"- 

We start with the iV-body Hamiltonian 



H 



1 ^ 



'-AiY + V'(fl,X2, ...,XAr) 



It turns out that we are not free to choose the rest of 
Rij arbitrarily. If we want to have a well-defined angular 
momentum decomposition in terms of the new variables 
{wi}, we will need to ensure that the resultant Hamilto- 
nian is of the same form as the las t line of Eq. |B1| This 
is because the kinetic term in (Bl I can be written as 



1 ^ 

— T 



-4/2 " 



1 I l2 /- <9 



(B2) 



with eigenstates labeled by angular momentum m and 
LL index n: 



|n, m) 



(6t) 



m+n 



(at) 



■\/(m + n)! "v/nT 



|0,0). 



(B3) 



This can be seen from how the angular momentum op- 
erator exist as part of the kinetic single-particle Hamil- 
tonian. With the second-quantized operators defined by 
(particle index i suppressed) 



1 / z „, a 



1 ^ 

— F 



' dxi 



c 2 



eB Xi 



+ V{xi,X2, ■■.,Xn) 
Q2 



2m 

+ V{xi,X2, ...,Xjv), 



(Bl) 



where Zi = Xi-iyt, Ai = B{-y,, x,,0)/2 and Ib = 
The symmetric gauge has been used so that the Hamil- 
tonian eigenstates will be conveniently labelled by angu- 
lar momentum. However, the results that follow will be 
gauge-invariant. For now, we will not make any assump- 
tion about the form of the interaction potential V. 

Haldane's original procedure was to first seperate this 
Hamiltonian into a CM part and relative part, and then 
project the relative part into different angular momen- 
tum sector^. The same will be done here, but with N 
particles instead of two. Define a change of coordinate 



= R 



with wi — {zi + Z2 + - ■ ■ + zn)/N, i.e. the CM coordinate. 
Any part of the resultant Hamiltonian depending only on 
wi will not affect our PP expansion. 



1 



2lf 



d 



V2\21b dz 
the angular momentum operator is 



L = h\ z— - z — 
oz oz 



To expand in terms of angular momentum eigenstates, 
the coordinate transform Ra must leave the form of each 

invariant, i.e. IziP must 



Bl 



term of the last line of Eq. 

transform into a sum of similar quadratic terms, etc. 

Denoting z = (zi, Z2, z^r)-^, dz = 
(sir' si;' sl;^)'^' likewise for the w's, the 

various terms transform as 



E 



dzidz 



dz'^dz ^ duF[RR^]dw, (B4) 



^ |z,|2 = ^ w'^[{R-^)T[R-^)]w = w^[RR^]-^w, 

(B5) 

E -^^T^- = z'^dz = w^[{R-^fR^]dw = w'^dw. (B6) 
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The Hamiltonian retains the same form if RR^ is diag- 
onal. If we regard i? as a rotation matrix, we see that 
this condition is satisfied whenever R maps an orthogo- 
nal basis to another orthogonal basis. Hence an allowed 
R consists of mutually orthogonal rows. As a simple ex- 
ample, the R matrix for 2 particles satisfies the condition 



1/2 1/2 



1 



-1 



according to the CM coordinate wi 
the relative coordinate 'W2 = Z2 — zi. 



{zi + Z2)/2 and 



2. The explicit form of the Hamiltonian 
transformed into total relative coordinates 



The next step is to explicitly find the coefficients of 
the transformed Hamiltonian. Since we are interested 
in a PP expansion in the total angular momentum, we 
define the total relative coordinate 



N-l 



W2 



Zn) = 



E 



N-l 



N-1 



- ZN. (B7) 



The other coordinates can be arbitrarily defined as long 
as they are orthogonal to W2 and wi = J2i' ^i- With 
this choice, the diagonal elements of RR^ are 



Ai = iV,A2 = 



N 



N-r 



Hence the kinetic part of the Hamiltonian becomes 

N r o2 



kin 



2m 

1 

2m 

+ {W2 



2 ^ 

^ * dwidwi 



d _ 
dwi 

d 
dwo 



+ 



1 



dW2dW2 4/2 



^1 



\W2\ 



W2T. j 

aw2 



(B8) 



where Irei = Ib\P^ = /sy is the effective "magnetic 

length" for the total relative coordinate. Only the terms 
corresponding to the total relative coordinate are shown 
in the second lino. In general, the diagonal elements of 
RR^ Ai, A2, A3, Ajv define a set of effective magnetic 
lengths li — Ib^TK- 



3. Derivation of the N-body pseudopotential 

We are now set up to find (mi, TOjv|l^|mi, ...,mN), 
the projection of an interaction potential 
V(wi,W2, ■■■,wpf) onto the angular momentimi sec- 
tor in the LLL. This projection is of course dependent 
on w = Rz. For the wi and W2 previously defined 
as the CM and total relative coordinates, mi and m2 
correspond to the CM angular momentum and total 
relative angular momentum, respectively. To evaluate 
this matrix element, we Fourier transform to shift the 
coordinate dependences onto an universal exponential 
factor: 



U 



= (47r) (mi, ...,mjv|F(u;i, ...,wjv)|mi, ...,mN) 



= (47r)^(mi,...,m^| j ^^g^ | y(fci, fc^) [] e^'^""' |mi, m^) 



= (4-)^ n 



{2^f 



\{{mj\V(ki,...,kN)e'''^-^^\mj) 



3 13 

''^ '^3 ^3 ) 



(B9) 
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The momentum-space potential V in the last line can be taken out of the expectation value since the momenta labeled 
by kjS are regarded as complex numbers. For each j, 



(TOj|e'''''"^|TO 



— ymj\e 2^ ' '\mj) = {jnj\e^'^ ' ^'e-^^ ^ \^j) 



s=0 



(s!)2(to — s) 



(BIO) 



The terms containing the aj and operators in the third line reduce to unity because the states are already defined 
to be in the LLL. Use has been made of the Baker-Campbell-Hausdorff formula in producing the factors of e^''^' 'j^"*. 

I 



The entire LLL pseudopotential reads 

jy-mim2...mjv 

= (mi, ...,mN\V{wi,W2, ...,WN)\mi, ■■■,mjy) 



N 

n 



V{ki, ...,kN), 



where m 



This reproduces the familiar result U^^{k) 
for N = 2. For TV = 3, [/"(fc) = (ffc^^l) 
is the total relative angular momentum characterized by 
W2 = ~ -^a) + (^2 — ^s))- The latter will be used 

extensively in the calculations of the next section. 



(Bll) 

where T^(fci, fc^r) is the Fourier transform of 
V{wi, ...,wm)- We focus on the W2 degree of freedom. 
Let I to) denote the state where TO2 = to and all other 
TOi — 0, i.e. the state with total relative angular momen- 
tum TO. Then the TOth PP for a translationally invariant 
interaction is 



N 

n 



4. Discussion and Generalizations 

We can determine how the effective magnetic length Ij 
should be rescaled without assuming the explicit form Wj 
defined in terms of the ZiS. By examining the diagonal 
elements of RR^ , we find that 



L — If 



\ 



N-l 



^ RijRi 



j<k 



xL,-, 



V2(iV-l) ;^(fc2,A:3,...,fc^). (B12) 

The integral over ki has been omitted since V does not 
depend on the CM coordinate wi. Also, the rest of the 
Laguerre polynomial factors have disappeared since Lq = 
1. We stih need to define f/"^-"" (/cs, fcg, Aiat) such 
that it corresponds to a PP 

(toi, ...,mN\V{wi,W2, ■.■,WN)\mi, ...^mpf) 

that is nonzero only at the simultaneous set of angu- 
lar momenta m2, ...,tojv. This is done by exploiting the 
orthonormality relation 2qe~'^ Ls{q'^)Lt{q^)dq — 5st- 
Switching to polar coordinates , we see that 




(B13) 



where Vq is again a constant with units of energy. If 
we want a PP that has no angular momentum on the 
spurious degrees of freedom , 1114 , ... , we find 



[/"(fc) = VoL,, 



eilN 

2{N -1) 



(B14) 



for i > 1. For i = 1, /i = VNIb- Here, wi is the CM 
position and li is the effective magnetic length for the 
CM angular momentum. 

(Bill can also be extended to cases beyond the LLL. 



There, the a and terms in the third line of Eq. BIO 



will not yield unity. If we consider the case where each 
particle occupies a specific Landau Level, we will have to 

first calculate expressions such as e V2^'^j V2^°'^ before 
making the change of coordinates from z io w. This is 
because the positions of the particles are indexed by z, 
not w. 

To begin with, we rearrange the exponential factor 
in the Fourier transform e*'^^''"^ = e*^'^^^^')^' so that 
it depends explicitly on the ZiS, albeit with modified 
k = (kjRji), summation implied. After some algebra, 
we find 



/ / 1 ka\ kaA \ \ 

(n|e^ ^e^^ '\n) 



ilsk 



where n and n' denote the initial and final Landau levels 
of the particle. When n = n' , we just have L„(fc^Z^/2). 
With the states being reduced to the LLL, we proceed 
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as in (Bll), arriving at the general formula for the 
PP between N particles initially at LLs ni,n2, ■ ■ ■ ,njsf 



mapping to LLs n'l ,n'2, ■ 



'ATI 



with angular momenta 



TOi, m2, . . . , m^r associated with the coordinates Wi 



n\ ...n'j 



..rriN 

;ni ...TIN 



= ("i, ■■■,n'j^;mi, ...,mN\V{wi,W2, wjv)!^!, t^at; Wi, ...,mjv) 



N 

n 



N 

n 



—il) 




\{Rukif\ 
2 



«JVJ 



(B15) 



with the effective magnetic lengths Ij as before. 

If V is translationally invariant, it does not depend on 
Wi — (zi + • • • + zi\[)/N, and the ki integral produces 
a delta function S{ki). Hence, as is usually the case, ki 
should be excluded from all sums in Eq. B15 As an il- 



lustration for N ~ 2 bodies with interaction independent 
of the CM, only the k2 integration survives. We have 



Riiki = ±k2 and V = V(k2)- mi, the CM angular mo- 
mentum, is irrelevant for the interaction, so Vmi — 5mi,o- 
We can also deduce this result from the orthogonality of 
the Laguerre polynomials: if we impose the further re- 
striction n.i = for all i, i.e. particles stay in their 
respective LLs, Eq. |B15| reduces to 



j d^ke-'^B^il^illkDhr,, 



l%k2 



V{k2 



(B16) 
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Appendix C: Derivation of 2-body Fermionic 
Pseudopotential Hamiltonians on a Cylinder 



Here, we present the details of the derivation of the 
two-body PPs C/™. From Eq. [21] a special case of 
Eq. B16 the potential C/™ that is nonzero only for 



the relative angular momentum sector m is given by 
U^ir-r') = VoL,ni-pBV^)S^{r-r'), where Vo is a con- 
stant with units of energy. We find its LLL Landau gauge 

by projecting onto the 

_ e '° ' 



basis matrix elements Uli 



basis wavefunctions ipn (t) 



where Ly is the circumference of the cylinder: 



Vr 



nin2n3n4 ^ 

= T / 15 / '^'^'^''''^™(9'^l)e^'-^''"'^'Vi,(rX(r')^„3(^')^n.(^) + antisymm. (CI) 



r 



is a dimensionless ratio that is small the form 



Recall that k — 
in the limit of large magnetic fields. The two types of MT 
symmetry constraints mentioned in Section V are man- 
ifest in the above expression as (i) the CM conservation 
which corresponds to rii + 112 = ^13 -I- ?t.4, and (ii) one- 
dimensional translation symmetry, which is the invari- 
ance of C/™ under — >■ -I- a, where a is an integer and 
i — 1,2,3,4. CM conservation must be present because 
the J dy and J dy' integrals produce delta functions of 



and 



2tt 



(714 -ni) + qy 



27r 



(713 - ^2) - qy 



The CM conservation condition ni + n2 — 723 + 714, i.e. 
n = n' , appears after we combine these two delta func- 
tions. Since the CM of each LLL wavefunction occurs 
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along X = nlsn oc n, we see that the CM of the parti- 
cles must indeed be equal before and after a two-body 
hopping. 



We explicitly resolve the MT symmetry via the trans- 
lation Hi ^ rii + l, as Fourier terms cancel upon ~^ 

■>Jji — ^i^y This is exactly a magnetic translation 
under which the system in a magnetic field B — -M- is 
expected to be invariant. To continue the calculation, 
reduce C/™ to the / (Pq integral 



[7™ - ^° ^ 
V° L, 



«^/| / d'q5{qy + =-{h-l2))L^{lW)e 

Ly 



2tt 



antisymm 
antisymm. (C2) 



The second line follows from the fact that qy is constrained to be — ^ 
disappears from the expression, as expected from MT symmetry. A closed form expression for J7'" is given by 



^"(^2-^1) = j^ik-h)- Note that n completely 



U, 



1 I P 

m! 

^ pUm -p)l ^ ^ 

p=0 ^ ^ ^' j=Q r=0 



T{p-j + r + l/2)2"+J-'-3(iK)20-'-i) 



(C3) 



The lengthy expression above can be factorized into the 
form U[^i^ = gK^h'llbll, g = 4yo(27r)3/2, where 



p=0 



(-2)3p^^(«02p^(2j + i)! 

{j-py.{2p + iy. 



(C4) 



This result can be proven by induction. The first few 
&2J+is are 

1 2 ,2 



6f = ^(-3 + 4K2;2);g-K^i^ 



6f = 77gy(15 - 40k^;2 ^_ iQuH^e-''^'^ 




FIG. 8. Graphs of b'^^+^ from Eq. [cijfor j = 0, 4, 10, 17 with 
K = |. As J increases, the main region of contribution of 
6^^^^ shifts in the direction of larger \l\. 



bJ = —=(-105 + A20k^P - 336kH^ + MK^fye''^''^'' 



=0. 

Certain 6s are depicted in Fig. [8j Note that the ?7™ 
operators can always be decomposed into the product of 
two 6™ operators that are m degree polynomials of li 



Appendix D: Orthogonality of [/" 



The FQH PPs U"^ are not systematically orthogonal 



and ^2- These polynomials have the physically relevant once we place them on a cylinder or torus. We will study 



property that (i) the of higher m are "localized" at 
larger values of I and (ii) that they are orthogonal in the 
limit of K — > before we enforce the L^Ly periodicity in 
the x-direction of the cylinder. This is further explained 
in Appendix [D] 



the physical origin of this fact and make quantitative es- 
timates in this appendix. While we are only concerned 
with quadratic cases of ~ Ly in this paper, it will be 
instructive to investigate how deviations from orthogo- 
nality depend on general and Ly. 
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1. From plane to cylinder 

When we compactify the plane into a cyhnder, the 
notion of relative angular momentum is no longer well- 
defined. In a fixed LL, the relative angular momentum 
is proportional to the interparticle distance which can 
only be meaningfully defined when the latter is much 



smaller than Ly. Recall that the effective k 



27rif 



■f- after mapping to the FCI system. We thus expect 
orthogonality to occur only in the limit of large Ly, or, 
equivalently, small k. To proceed, we evaluate the overlap 
elements of two PPs VJyj and V^"), and show that the 
latter are orthogonal for sufficiently small k. 



hi 



cyl (i,;2 cyl li,l2 

J2j^'lJ dh'Siqy + ^{h-h))Siq[^ 
hh ^ 

E / dly 



h-h 



dqx 



dqx I dq'^S{qy ^ 
27r 



qy)Lm{lBq )Ln{lB<l' )e 



Ly 



dq'Jiqy 
2tt 



ih - l2))LUlW)Ln{lWy'''-''' 5{qx + q'x) 



\h — h 

TT J 2{qlB)d{qlB) (^Y. Hiqyh) + Kih - h))^ L„,{lW)Ln{lW)e''B'f 



(Dl) 



Two approximations have been made above. From 
the third last to second last line, we replace the non- 
rotationally invariant integral over J dq^dqy by the rota- 
tionally invariant integral tt J 2qdq. From the second last 
to the last line, the delta function sum in the large paren- 
theses was replaced by unity, i.e. taking the limit where 
the range of li — h (and hence q) tends to infinity. These 
approximations become exact when the discrete qylB be- 
comes a continuum, which occurs precisely when k — 0. 
Indeed, this agrees with the physical intuition that the 
relative angular momentum becomes a well-defined quan- 
tity when Ly is large. From the viewpoint of polynomial 
orthogonality, we see that the orthogonality of the ?7™j is 
respected as much as as the integral over Laguerre poly- 
nomials is allowed to be made continuous. Specifically, 
the orthogonality of the Laguerre polynomials is exact 
only for the continuum q^ = q^ + Qy and not for discrete 
points specified by Al — li— 12 whose separations do not 
vanish unless k = 0. 



r 



space. Nonorthogonality is expected when there is sig- 
nificant overlap between these images. Let us obtain a 
bound by finding the conditions where |f/^j|, whose ex- 



plicit form is given by (C4), is not negligible, with li,l2 



being of the order of L^Ly. At such values, 



cyll\l2 



'2w 



\2m 



f{w) 



(D2) 



where w = nL^Ly. The function f{w) exhibits a rapid 
Gaussian decay beyond w ~ 2y'rn. Hence we expect the 
PP [7™ to be nonorthogonal when 



K < 



LxLy 



1. FVom cylinder to torus 



The compactification of the cylinder into a torus intro- 
duces a periodicity in the La;-direction. This introduces 
periodic copies of f7™j in {7^™ , each displaced from an- 
other by {L^Ly, LxLy) or {L^Ly, -L^Ly) sites in - I2 



which implies that 



Lx < 2\fm 



(D3) 



a bound well- verified by calculations in the following sub- 
section. 
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FIG. 9. (Color online) The eigenvalues of the overlap matrix 
of ,U'^ jU"^ and as a function of the system size 
L — — Ly. Indeed, we observe orthogonality when L > 6, 



in excellent agreement with the bound 2\/9 — 6 from Eq. D3 



Numerical results for the orthonormality of the 
fermionic pseudopotentials 



The orthonormality of the U™s can be studied quan- 
titatively through their overlap matrix defined by 



according to Eq. 36 When the [/™s are orthonormal, it 
holds Mmn = id. If the [/™s are not orthonormal while 
they span an orthonormal basis, the eigenvalues of M^n 
will still be unity since we can find a unitary transforma- 
tion where Mm„ is diagonal. The U"^s are overcomplete, 
however, the spectrum of M broadens and yields eigen- 
values deviating from this limit. In Fig. [9] we plot the 
eigenvalues of the overlap matrix for the first few PPs as 
a function of L. Orthogonality is hence broken when the 
eigenvalues differ from unity. Indeed, we observe that or- 
thogonality improves with system size, in agreement with 
the conclusions of the preceding subsections. 



Appendix E: Three-body bosonic Pseudopotential 
calculations in the LLL Landau gauge basis 

Let U"^ be the bosonic 3-body CM-invariant potential 
which is nonzero only in the sector of total relative angu- 
lar momentum m. Such 3-body U"^s will be useful as a 
basis through which arbitrary CM conserving potentials 
can be expanded. In this appendix, we will show the de- 
tailed derivation of C/™ in the basis of LLL Landau gauge 



eigenf unctions, starting from Eq. CI 



M„ 



{U"',U") 



According to Eq. B14 we have to replace the magnetic 
length l"^ by |?^. This, however, causes no difference 
when m = 0, since Lq = 1, so the results in Refs. [351 
and [55] for [/° could have been obtained without this 
modification along our GHPs. 

After performing the J (Pr^ integrations according to 
(D4) Eq. |CH we arrive at 



where 



-j^-^I{q,ni,ne)I{p - q,n2,n5)I{-p,n3,n4)L,n i^-{p-qfllj , (El) 

I{q,n, n') = 2Tr6 (^qy + ^ e-'««'/4g»Ki«(n+n')9./2g-K^(„-«')V4_ 

Here, T,a = JW ^permimn^ns) sgn(perm)^ Eperm(»6n5n4) sgn(perm)-^ where / 1 for fcrmions and / = for bosons. 
perm refers to the permutations of the n^s while sgn{perm) = ±1 depending on whether the permutation is even or 
odd. Up to now, the expression is easily generalizable to any number of fermions or bosons by increasing the number 
of integrals and I{q,n,n')s. If we further restrict ourselves to the 3-body bosonic case, we can simplify it to 



"-^nin2n3n4"5"6 

g-«;^([("i-«6)^ + (n2-"5)^-l-(n3-«4)^]/4 



^ f d q^d p^^ ^^ /3^^ qfll^j e~'s(p''+9'+(p-9)')/4g«isA'i9/2g«/BJV2P/2 



(x5]P„(ni,n2,n3,n4,n5,n6;«2)^e-^(^?+^^^^+^^) 



(E2) 



where Ni = rti + rtg — n2 — ns, = ^2 + — — 714. P,„(ni, 712, "-s, "■4, "-s, n^; k?) is a potentially complicated 
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polynomial whose form depends on m. We define the 
conserved CM "Landau level wavefunction index" R via 
3R mod N = rti + 71,2 + = ng + ri5 + 714. The existence 
of i? is a conseqence of the total CM conservation of U™. 



1. m=0 for bosons on a torus 

In this case, Lq = 1, and PQ{ni, n2, n^, ^^4, ^^5, Uq; k^) = 
1. Since each of the n^'s is also defined modulo N (but 
constrained to sum to 3i? as shown above), we make the 
replacement 

g-4^((i?,-ni)^ + (fl-n2)^ + (-R-n3)^) _^ 
s,t 

(E3) 

and likewise for the identical factor involving 77.4,77.5,77.6. 
Hence the C/'" factorizes into a product of nonlocal op- 
erators 



Hence we have 



where 



E 

X]; ni=3Ra\odN 



L-ni *~^n2 '-'^^3 



E 

^^ 1 + TT. 2 + 3 — 3 i? m o d 



E' 



L-ni *--n2 1 



(E4) 



with Wst = Z]i"-'i - J2i<jKn'ji "i = 'T-i + s^, n2 = 
77-2 + tA^ and 713 = 773 — N{s + t). This result is identical 
to that in Ref. [55l 



2. m=l for bosons on a torus 



C/i = 0. 



(E6) 



This agrees with the result from Ref. 69, in that there 
is no PP of total relative angular momentum 777 = 1 for 
bosons. This is because the CM invariance of the inter- 
action precludes any symmetric wave function of total 
degree 1. 

3. m=2 for bosons on a torus 

After some algebra, 

P2 = ^(1 - 3^2(772 - Rf)il - 3Ac2(n5 - Rf). 



Since 



^(772 - = = ^(772 - 2n2R + 

(T <7 

= \j2nl~2R^ + R^ 



i i<j 



we have 



where 
bR = 



E 

Y,ini=3RmodN . 



2^2 



■-•ni L-n2 



(E8) 

with Wst = J^i^'i - Y.tKjK'^'j^ "-i ^ ni + sN, n'^ = 
772 + tJ^, and 773 — 773 — N{s + 1) as before. 



4. m=3 for bosons on a torus 

This is zero for 2-bosons, but not for 3-bosonJ^. In- 
deed, 



In this case, the 1st Laguerre polynomial gives a factor 

4\ 



1 — |(k2(772 — 775)2 +Px), and Yla evaluates to 



E ^3 = E - Rf{n, - Rf = n(77, - R). 



^ Pi - ^ 3K2(n2 - i?)(775 - i?) 



3k^ 




-R 



-R 
(E5) 



Hence 



where 



cx ^ 5|j6 
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6r = 



E 



■n. 1 + n 2 + n 3 — 3 i? m o d A'' 



^(ni + sA^ - i?)(7i2 + W - R){n'i - (s + ^)iV - i?)e" 



^rii ^712 '-"H-a ■ 



(E9) 



5. m=4 for bosons on a torus 

After the smoke clears, we find 



with Wst = Y^i^'l - Y.tKjK'^'j^ = ni + sN, n'^ = 
n2 + tN, and n'^ — ~ N{s + 1) as before. 



R 



where 
bR = 

E 

ni=3RniodN 



Cm ^712 '^ns 1 
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